Spin polarized electron transport near the Si/Si02 interface 
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Using long-distance lateral devices, spin transport near the interface of Si and its native oxide 
(Si02) is studied by spin- valve measurements in an in-plane magnetic field and spin precession 
measurements in a perpendicular magnetic field at 60K. As electrons are attracted to the interface 
by an electrostatic gate, we observe shorter average spin transit times and an increase in spin 
coherence, despite a reduction in total spin polarization. This behavior, which is in contrast to the 
expected exponential depolarization seen in bulk transport devices, is explained using a transform 
method to recover the empirical spin current transit-time distribution and a simple two-stage drift- 
diffusion model. We identify strong interface-induced spin depolarization (reducing the spin lifetime 
by over two orders of magnitude from its bulk transport value) as the consistent cause of these 
phenomena. 

PACS numbers: 85.75.-d, 72.25.Dc, 72.25.Hg, 85.30.Tv. 



Recent advances in the development of techniques for 
electrical injection and detection of spin-polarized elec- 
trons have revolutionized semiconductor spintronics.Q 
Among the materials studied so far, silicon (Si) stands 
out due to its weak spin orbit and hyperfine effects, 0, Q 
effectively insulating electron spin from relaxation mech- 
anisms and leading to long spin lifetimes (e.g. > 500ns 
at 60K) in bulk.0] Although this long lifetime was sub- 
sequently exploited to demonstrate long distance spin 
transport of over 2mm in quasi-lateral devices,!^ the 
transport mode was still largely bulk in character, where 
the inversion symmetry of the diamond lattice potential 
protects spin degeneracy of the conduction band and sup- 
presses spin-orbit effects. 

At interfaces where bulk transitions to 2-dimcnsional 
transport, however, the lattice inversion symmetry is 
broken. This opens the possibility of the emergence of 
(weak) spin orbit effects including Bychkov-Rashba spin 
manipulation. [S| Besides the bare surface of Si, which is 
prone to contamination and is difficult to electrostati- 
cally gate,0| the buried Si/Si02 interface is of interest in 
exploring this physics experimentally because of its high 
conduction band offset with the insulating oxide and its 
technological importance in charge-based Si field-effect 
devices. 

Although it is widely expected that the approach 
to two-dimensional transport tends to preserve spin 
polarization as a result of limiting orbital degrees of 
freedom,^ here we show (using gate-tunable lateral Si 
spin transport devices and a corresponding transport 
theory) that attracting bulk conduction electrons to the 
thermally-grown Si/Si02 interface causes a massive spin 
lifetime reduction. These results are important for the 
development of gated semiconductor spin transport de- 
vices, such as the canonical Datta-Das transistor, Q and 
may have a similar physical origin to effects seen in spin 
transport in graphene/Si02 devices, where the observed 
spin lifetime is many orders of magnitude lower than 
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FIG. 1: (a) Schematic view of the electrostatically-gated lat- 
eral Si spin transport devices, showing electrical connection 
configuration, (b) Plan-view micrograph of a typical device. 



expected. 

Our lateral spin transport devices utilize ballistic hot- 
electron spin injection and detection techniques, simi- 
lar to previously-reported vertical-transport devices. 0, 
H, m, [13, Q, [H, Ell However, the lateral geom- 
etry used here required modified fabrication methods. 
First, to make a semiconductor- metal-semiconductor 
spin detector, ultra-high vacuum metal film wafer bond- 
ing was performed with two different silicon-on insula- 
tor (SOI) wafers using a thermally evaporated NigoFe2o 
(4nm)/Cu (4nm) bilayer; one SOI wafer consists of 
3/im single-crystal (100) n-type silicon for the detec- 
tor hot electron collector and the other of 10/im single- 
crystal (100) undoped sihcon for the transport chan- 
nel. The rather large thickness of lO/im was chosen 
to achieve wide proximity control of spin-polarized con- 
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duction electrons to the Si/Si02 interface below it us- 
ing the p-type Si handle substrate wafer as an electro- 
static gate. Conventional wet-etching methods, along 
with photolithography, were used to pattern the spin 
detector and the subsequently exposed undoped sili- 
con layer was Ar+ ion-milled to eliminate a conduc- 
tive silicide on the surface, which could be detrimen- 
tal to spin transport. [3| On the exposed undoped sil- 
icon layer, a hot-electron spin injector was assembled 
with a ferromagnetic cathode tunnel junction having 
a structure of Al (40nm)/Co84Fei6 (10nm)/Al2O3/Al 
(5nm)/Cu (5nm)/undoped-Si (10 fim). A schematic di- 
agram of the completed lateral device structure and a 
plan- view micrograph are shown in Figs. [Ha) and (b), 
respectively. The distance from the edge of the base (the 
Al/Cu metal layer under AI2O3) to the edge of the detec- 
tor is approximately 150 ^m, but the edge of the emitter 
tunnel junction is receded from the base edge by approx- 
imately 20/im. XRD measurements confirmed that elec- 
tron spin transport from the injector to detector is along 
the [Oil] crystalline axis. 

To perform electron spin transport measurements, a 
tunnel junction emitter voltage V^; larger than the Cu/i- 
Si injector Schottky barrier height is applied to inject 
spin-polarized hot electrons ballistically into the 10 /im- 
thick undoped silicon channel. Because the lateral elec- 
tric field caused by the accelerating voltage drop Vci is 
screened by the equipotential injector base layer, these 
electrons first must diffuse to the channel region and are 
then carried primarily by drift toward the detector. Af- 
ter coupling to hot electron states in the Ni8oFe2o, the 
ballistic component of the charge current is filtered by 
the asymmetry between spin-up and spin-down mean- 
free-paths in the ferromagnet, so that the number of 
electrons coupling to conduction band states in the n- 
type Si detection collector (comprising the signal cur- 
rent Ic2) will indicate the final spin polarization of the 
transported electrons. Low operating temperatures are 
required to freeze out unpolarized thermionic currents 
over the Schottky barriers, which would otherwise dilute 
the small ballistically-injected spin polarized current. 

Two types of measurements were performed using ex- 
ternal magnetic field sweeps, first in an in-plane geometry 
and second in a perpendicular field geometry at 60K with 
constant Ve =-1.3V. The resulting spin-valve and spin- 
precession data at constant Vc'i=20V are shown in Fig. 
[2] for gate voltage bias Vg=-4V, where injected electrons 
are pushed away from the Si/Si02 interface, ((a) and (b), 
respectively), and Vg=-I-3.5V, where they are attracted 
to the interface, ((c) and (d), respectively). 

A comparison of the spin-valve data in Fig. [H^a) and 
(c) reveals similar device response: As the external mag- 
netic field sweeps through the in-plane coercive field of 
the magnetically softer Ni8oFe2o, an antiparallel injector- 
detector magnetization configuration is obtained, and the 
signal current determined by projection of spin orienta- 
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FIG. 2: (a) Spin- valve measurement in an in-plane magnetic 
field and (b) spin-precession (Hanle) measurement in a per- 
pendicular magnetic field under negative gate biasing condi- 
tions which push transported electrons away from the Si/Si02 
interface, (c) Spin valve measurement and (d) spin precession 
measurement under positive gate biasing conditions which 
pull transported electrons toward the Si/Si02 interface. All 
measurements use emitter voltage Ve =-1.3V and accelerat- 
ing voltage Vci =20V at 60K. Arrows indicate direction of 
magnetic field sweep for the appropriately colored data. 



tion on detector magnetization Ic2 is reduced. This is 
due to a decrease in ballistic mean-free-path in the detec- 
tor ferromagnet, and is identical to the behavior of ver- 
tical spin trans por t devices utilizing the same detection 
technique. (3. [Bl. Hi. T^ . isl 14, 1^ However, marked 
qualitative differences are clear from a comparison of the 
spin precession data in a perpendicular magnetic field, 
shown in Figs. E^b) and (d): although the magnitude of 
spin precession oscillation is reduced as the gate voltage is 
varied from negative in Fig. [U^b) to positive in Fig. [IKd), 
the coherence (e.g. the number of clear oscillations) is 
dramatically increased. Furthermore, the magnetic field 
period of these coherent oscillations is larger in Fig. [2td) . 

Since these oscillations are due to spin precession at 
angular frequency lu — g/jEB/h (where g is the electron 
spin g-f actor, /is is the Bohr magneton, B is the mag- 
netic field, and h is the reduced Planck constant) dur- 
ing transport time t, the period of the oscillations B2Tn 
caused by a change in the average spin orientation, can 
be used to recover t = 2tt?i/ gfiBB2Tr-^i^ Therefore, the 
larger oscillation period seen in Fig. [U^d) indicates a 
shorter transport time, despite the identical voltage drop 
and drift field between injector and detector. 

This behavior is in contrast to vertical Si spin trans- 
port devices, where shorter spin transport time leads to 
an increase in observed spin polarization due to a sim- 
ple exponential spin decay mechanism. 4] However, the 
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geometry here is very different; for example, the lateral 
extent of the injector (« 400/xm) and detector (« SO/xto) 
is expected to cause an unavoidable systematic source of 
decoherence (or spin "dephasing" ) because of a transit- 
length uncertainty leading to transit time uncertainty At 
and suppression of oscillations due to spin precession an- 
gle uncertainty A0 = ujAt. 17[ In light of this expected 
geometrical dephasing, the many coherent oscillations 
seen in Fig. [2jd) is surprising. 

One simple explanation for this observed behavior 
could be a spin lifetime on the order of the transit time 
t. This would suppress contributions to the magnetic 
field-dependent spin signal by electrons which arrive at 
the detector after long transit times, shortening the ob- 
served average transit time and the overall polarization, 
and decrease the width of the distribution and therefore 
the dephasing. However, previous studiest4|| have estab- 
lished spin lifetimes in bulk undoped Si using vertical 
transport devices at this temperature of at least 520 ns, 
and the apparent transit time here is « 2ns, so a drastic 
reduction in spin lifetime induced by the presence of the 
interface would be necessary to support this explanation. 

Fortunately, we have empirical access to the spin cur- 
rent transit time distribution through the spin precession 
data. This signal is a weighted sum of cos ujt contribu- 
tions (due to precession at angular frequency uj) from 
the projection of spin orientation onto the magnetization 
axis of the detector thin film.[i,[il,[i3,lil The weighting 
function P{t) is just the spin current transit time distri- 
bution at the detector so that our expected signal is given 
by: 



Ic2 oc j P{t) cosojtdt. 



(1) 



Because P{t < 0) = 0, we can extend the lower bound 
of integration to — oo in Eqn. [T] without consequence, 
and recover an expression equivalent to the real part of 
the Fourier transform of the spin current transit time 
distribution: 



Ic2 oc Re 



(2) 



Since the spin precession signal Ic2 is real, the transform 
has even symmetry in w, so the real part of the Fourier 
transform of the spin precession data will then recover 
the spin current transit time distribution. 

When we perform this transformation, spin current 
transit time distributions such as those for Vci =10V 
shown in Fig. [3l^a) are extracted. At a gate voltage 
bias of Vg'=-4V, the distribution is wide and peaked near 
10ns. However, as the gate voltage is made more positive 
and electrons are attracted toward the Si/Si02 interface, 
the resulting spin transit time distribution contracts to 
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FIG. 3: (a) Spin current transit-time probability distributions 
for different gate voltages Va as labeled, derived from Fourier 
transforms of spin precession measurements using Eqn. [2] for 
Vci =10V. Distributions are normalized by the total spin sig- 
nal Ic2 in a magnetically parallel injector and detector con- 
figuration, (b) Spin transit-time distribution functions calcu- 
lated using the two-regime drift-diffusion model described in 
the text for different spin lifetimes as labeled. 



lower transit times but always remains within the wider 
distribution. This behavior clearly suggests a suppres- 
sion of contributions from long transit-time electrons by 
fast depolarization at the Si/Si02 interface. 

To more quantitatively test this hypothesis, we simu- 
late spin transit time distributions using a 1-dimensional 
model where injected spin-polarized electrons (screened 
from the lateral drift electric field by the injector base 
layer) diffuse at the Si/Si02 interface until they reach 
the edge of the base, where they then are carried by drift 
and diffusion to the detector. Their total transit-time 
distribution is then given by the convolution of the tran- 
sit time distribution of the diffusion-dominated path un- 
der the injector base (where spins are injected a distance 
w < x' < w + W from the edge of the emitter base, 
where w is the exposed region of injector base under the 
emitter (here, 20 /im), and W is the length of the base 
covered by the emitter (here, 400 /im)) with the drift- 
dominated transit time distribution for transport in the 
L — 150/im-long channel between the injector and de- 
tector. The former (diffusion-only) distribution is given 
by a sum of all contributions (each the impulse response 
or "Green's function" of the spin diffusion equation[19]) 
from along the wide injector emitter to the beginning of 
the drift-dominated region: 
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where D is the diffusion coefficient and r is the spin life- 
time. The latter (drift-diffusion) distribution is just the 
Green's function of the spin drift-diffusion equation with 
absorbing boundary conditions at the detector x = L:(20j 
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FIG. 4: Total spin polarization, most probable transit time, 
and dephasing (given by the half- width half- maximum) , of the 
empirical spin transit time distributions ((a)-(c), respectively) 
and the simulated distributions ((d)-(f), respectively). 
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where v is the drift velocity. Neglecting any potential 
spin blockade effects '21'], the spin current distribution is 
then proportional to the diffusion current at the absorb- 
ing detector boundary, P{t) oc -D-£{pi (E)p2)\x=L- 

Simulated spin transit time distributions (for several 
spin lifetimes r =:2ns, 6ns, 14ns, and 22ns) are shown 
in Fig. [31[b) for comparison with the empirical distribu- 
tions in Fig. El^a) . The transport parameters (determined 
from fits of spin precession measurements from 350um 
undoped-Si vertical spin-transport devices at 60K)[12| 
are mobility ^ =5833 cm^/Vs and D =1500 cm^/s. 
We do not intend an exact fit of the empirical distri- 
butions; rather, we show that this model captures the 
salient qualities such as the similar peak-probability ar- 
rival time shift, polarization reduction, and transit-time 
uncertainty (At) suppression with spin lifetime reduc- 
tion. The details of the empirical distributions are likely 
difficult to fit with a 1-dimensional model like the one 
we have employed here, since the electrostatics in the Si 
transport channel will likely couple the otherwise inde- 
pendent voltage bias variables Vci and Vg- 



Further support for the conclusion that the interface 
contributes to a reduction in spin lifetime is provided by 
a comparison of empirical and simulated spin transport 
distribution characteristics shown in Figs. IHa)-(c) and 
(d)-(f), respectively. As the drift electric field increases 
with increasing accelerating voltage bias Vci , the empir- 
ical spin polarization saturates, transit time decreases, 
and spin dephasing (here, the distribution half-width at 
half maximum) decreases. At any constant Vci, increas- 
ing Vg tends to decrease all these parameters. Analysis 
of the simulated distributions shows a very similar be- 
havior as shown in Figs. IHd)-(f). Remaining differences 
in the value of spin polarization between Figs. lUJa) and 
(d) can be accounted for by a non-ideal spin polarization 
at the injection point, and the large-Vcr separation be- 
tween transit-time uncertainty (dephasing) seen in Fig. 
mf) that is not seen in the experimental results in Fig. 
IHc) is likely due to the details of channel electrostatics 
as mentioned previously. 

The microscopic cause of the observed severe spin 
lifetime suppression close to the interface is not clear. 
Interface roughness and presence of charged impurities 
are expected to enhance momentum scattering and in- 
crease spin depolarization proportionately via Elliott 
scattering, [2^ but the experimental spin transit distri- 
butions do not show a shift of the low-transit-time edge 
to longer times that would result from such a decrease 
in charge-transport mobility. However, the random lo- 
cations of these charged impurities could cause a fluc- 
tuating Bychkov-Rashba field leading to an enhance- 
ment of spin relaxation. [2^ The possibility also remains 
that the broken lattice inversion symmetry at the inter- 
face causes a D'yakonov-Perel'-type spin depolarization 
mechanism, [2^] and our future studies using these de- 
vices will explore whether this is the case, or if differ- 
ent effects such as localized spin exchange with param- 
agnetic defects in the oxidef25] or dangling bonds at the 
interface[2^ are to blame. Regardless, this effect could be 
exploited for novel devic e op eration requiring electronic 
control over spin lifetime, [27] and could at least partially 
be to blame for an absence of evidence for spin injection 
in prior attempts at spin transport in Si. 
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